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ABSTRACT 

We estimate the two-point correlation function of dark matter haloes, with 
masses > 10^^ Mq, that have or not significant substructure. The haloes are 
identified with a friends of friends algorithm in a large ACDM simulation at two 
redshift snapshots {z = 0.0 and 1.0), while halo substructure is determined using 
an observationally driven method. We find in both epochs a clear and significant 
signal by which haloes with substructure are more clustered than those with no- 
substructure. This is true for all the considered halo mass ranges, although for 
the highest halo masses the signal is noisy and present only out to ~ 20h^^ Mpc. 
There is also a smooth increase of the halo correlation length with increasing 
amplitude of the halo substructure. We also find that substructured haloes are 
typically located in high-density large-scale environments, while the opposite is 
true for non-substructured haloes. If the haloes found in high-density regions have 
a relatively earlier formation time, as suggested by recent works, then they do 
indeed have more time to cluster than haloes, of a similar mass, which form later 
in the low-density regions. In such a case one would have naively expected that 
the former (earlier formed) haloes would typically be dynamically more relaxed 
than the latter (later formed). However, the higher merging and interaction rate, 
expected in high-density regions, could disrupt their relatively relaxed dynamical 
state and thus be the cause for the higher fraction of haloes with substructure 
found in such regions. 

Subject headings: cosmology: theory - dark matter - galaxies: haloes - galaxies: 
formation - large-scale structure of universe - methods: N-body simulations. 
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1. Introduction 

In the current cold dark matter picture, structures form via gravitational amplifications 
of initial density fluctuations. The dark matter (DM) haloes form hierarchically by the 
aggregation, via gravitational interactions, of small collapsed structures that merge to form 
larger ones, eventually forming clusters of galaxies. 

Observationally, galaxy properties show a significant variation depending on their local 
and large-scale environment (eg. Dressier 1980; Gomez et al. 2003; Boselli & Gavazzi 2006). 
Since galaxies form within DM haloes, the "parent" halo properties could influence or even 
determine the galaxy properties. Therefore, many theoretical studies have investigated in 
detail the inter-relation of the DM halo properties, among which the halo formation time, 
their assembly histories, their central concentration, their clustering and their location in 
the cosmic web. The first studies (eg. Lemson & Kauffmann 1999) found no significant 
environmental dependence of the various DM halo properties. However, a re-interpretation 
of their results and a wealth of new studies point in the opposite direction. 

Halo formation time has been found to correlate with halo concentration and with the 
number of sub-haloes (eg. Navarro, Frenk & White 1996; Jing 2000; Bullock et al. 2001; 
Wechsler et al. 2002; Zhao et al. 2003; Zentner et al. 2005). Sheth & Tormen (2004) 
showed that haloes in dense regions form at slightly earlier times than similar mass haloes 
in lower density regions, a result confirmed by Zhu et al. (2006) and Harker et al. (2006) 
for massive haloes as well. Gao et al. (2004) and Gao & White (2007) extended the halo 
"assembly bias" investigation by considering different parent halo properties, among which 
concentration, spin and subhaloes mass fraction. 

Furthermore, the amplitude of the DM halo two-point correlation function has been 
found to depend strongly on the halo formation time. Relatively low-mass haloes (^ lO^^/i^^M©), 
that assembled at high redshifts, are more clustered than those that assembled more recently 
(Gao, Springel & White 2005; Wechsler et al. 2006; Wang et al. 2007; see also Berlind et al. 
2006), an effect which strengthens with decreasing halo mass and with decreasing halo sepa- 
rations. However, Wetzel et al. (2006) and Jing, Suto & Mo (2007) found that very massive 
haloes, which have formed recently, appear to cluster more strongly than older haloes of the 
same mass, which is the opposite than what has been found for lower mass haloes (see also 
Wechsler et al. 2006). 

Regarding the dependence of halo clustering on subhalo occupation number, Wechsler et 
al. (2006), studying only large host haloes, found close pairs of haloes having above-average 
number of satellites. Gao & White (2007), using the Millennium simulation (Springel et 
al. 2005), demonstrated that haloes with a relatively large mass fraction in sub-haloes are 
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more strongly clustered, confirming Wechsler et al. (2006). Thus halo occupation by galaxies 
could be a function of environment. This appears to be in agreement with some observational 
cluster studies, to the extent to which halo-occupation is related to the halo dynamical state. 
For example, dynamically active clusters (having significant substructure) have been found 
to be preferentially located in dense large-scale environments (Schuecker et al. 2001; Plionis 
& Basilakos 2002), while they are also more clustered with respect to clusters with weak or 
no substructure (Phonis & Basilakos 2002). 

Further on the observational side, Blanton et al. (2006) found a slight dependence 
between galaxy colours and measures of the surrounding environment on various scales. 
Berlind et al. (2006), searching for a dependence of the SDSS galaxy groups clustering on a 
second parameter, independent of mass, found a correlation between the large scale bias of 
massive groups and their central galaxy (g — r) colour as well as hints of a possible connection 
between the later and DM halo age. Also Croton, Gao & White (2007) studying the relation 
between the clustering of DM haloes and their assembly history, using a galaxy formation 
model applied in the Millennium simulation, found that the colour of the central galaxy in 
a DM halo, of a given mass, depends on the halo's environment. 

All these results put into question the simplest form of the excursion-set formalism for 
galaxy clustering (Bond et al. 1991; Lacey & Cole 1993), which predicts that the properties 
of galaxies, forming within the DM haloes, arc functions only of the DM halo mass, and 
not of other properties like its environment. Nevertheless, more elaborate schemes of the 
excursion set formalism have been recently proposed that allow more faithful descriptions 
(Sandvik et all. 2007, Zcntner 2007). 

In this letter we address a related issue; ie., the relation between DM halo clustering, the 
halo dynamical state and its large-scale environment. Note that we do not use the common 
approach of counting bound sub-haloes, within parent haloes, as an indicator of substructure, 
but rather an observationally driven estimator of the halo dynamical state which is based 
on the dynamical measure of halo substructure according to Dressier & Shectman (1988). 
Once the halo dynamical state is determined we compute the two-point correlation function 
for haloes with and without substructure and we also characterize their local environment. 



2. Numerical Simulation & Methodology 

The numerical simulation used in this work was performed using the GADGET2 code 
(Springel 2005) with dark matter only. We use a ACDM cosmology with = 0.3, Qa — 0.7, 
as — 0.9, h — 0.72, where and are the present day matter and vacuum energy densities 
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in units of the critical density, as is the present hnear rms amphtude of mass fluctuation 

in spheres of 8 Mpc and h is the Hubble parameter in units of 100km s~^Mpc~^. The 
simulation was run in a cube of size L = 500 Mpc, using 512^ particles. The particle 
mass is ~ 7.7 x 10^° Mq and the force softening length is e = 100 kpc. 

The DM haloes were identified using a friends-of-friends (FOF) algorithm with a hnking 
length I — 0.17 times the mean intcr-particle separation. Given the purpose of this work, 
we only use haloes with at least 130 particles, ie., with masses greater than IO^'^^Mq, which 
results in a sample of ~ 58000 haloes at 2; = and ~ 32400 at z = 1. Note that the halo 
finding algorithm used provides unique haloes that are not sub-haloes of any other "parent" 
halo. 

Since we wish to investigate the possible correlation between clustering and the dynam- 
ical state of DM haloes, we use the Dressier & Shectman (1988) algorithm to estimate the 
amount of halo substructure. Details for a recent application of this method to simulation 
data can be found in Ragone-Figueroa & Plionis (2007). Briefly, this method determines the 
mean local velocity, (vjoc), and the local velocity dispersion, uioc, of the nearest n neighbours 
from each halo particle i and compares them with the mean velocity, (V) , and the velocity 
dispersion, a, of the whole halo of N particles, defining the following measure: 

71 

Si = -[((vioc) - (V))^ + (aioc - af] . (1) 

Then the quantification of the halo substructure is given by the so-called A- deviation, which 
is the sum of the individual 6iS over all halo particles A^: A = J2n ^i/^- The larger the 
A-deviation the stronger is the halo substructure. This statistic depends on the number of 
nearest neighbours n which is used in the analysis (eg., Knebe & Miiller 1999) and on the 
number of particles used to resolve a halo as well. Such resolution effects were studied in 
detail in Ragone-Figueroa & Plionis (2007), who found a monotonic increase of (A) with the 
halo mass which is clearly due to resolution effects. However, they also found that within 
relatively small halo mass intervals, the corresponding sorted A-deviation distribution can be 
used effectively to separate substructured from non-substructured haloes, by dividing them 
in those having A-deviation above and below the corresponding median or some quantile of 
the A-deviation distribution within each halo mass interval. 

Therefore, the total halo sample was divided in 3 mass ranges and the analysis was 
performed in each individual subsample and for each of the two redshift snapshots {z — 
and 1.0). These halo mass ranges are: 

(a) 10^3 H-^Mq < M < 3 X lO^^h-^MQ, 

(b) 3 X 10^3 h-^MQ <M < 10^^ h-^MQ, 
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(c) M > 10^^ H-^Mq. 

We compute the DM halo two-point correlation function, ^(r), separately for the substruc- 
tured haloes (dynamically active) and non-substructured haloes (virialized) , within each 
mass range, and to emphasize their possible clustering differences we separate haloes having 
A-deviation values larger than the 67% and lower than the 33% quantile of the A-deviation 
frequency distribution function. 

The measured halo two-point correlation function is also fitted to a power- law: ^(r) = 
{r/roy, in the range 4^ 30 h"^ Mpc, using a minimization procedure. Further- 
more, we estimate the halo mass one-point over-density distribution function, dividing the 
simulation box in grid-cells of 5 Mpc size and compare the total DM mass (using 
haloes with M > lO^^h'^MQ) within each cell, Mj, to its mean value, M, ie., estimat- 
ing: SM/M = [Mi — M) / M . We then assign to each halo the 5M/M value that corresponds 
to the grid-cell in which it is located and derive the corresponding frequency distribution, 
separately for haloes that have or not significant substructure. We have verified that this lo- 
cal density estimator is equivalent with the traditionally used, estimated in spheres centered 
on each individual halo. 



3. Results and Discussion 

In Figure 1 we show the ratio, -R(r), of the two-point correlation functions of haloes 
with and without substructure, for the three mass intervals and for both redshifts epochs. 
The errors shown are the propagated quasi- Poissonian uncertainties, estimated by: 8^{r) ~ 
A/[l + C(0]/DD(r), where DD(r) are the halo-halo pairs within separation r ± 5r. It is 
evident that R{r) is significantly less than 1, indicating that the former haloes indeed have 
a significantly higher correlation function than the latter ones. For 2; = we have R{r)^ 0.8 
for all r's, while for z = 1 and for the most massive haloes this is true for separations ^ 20 

Mpc. Note, however, that the size of our simulation does not allow us to probe effectively 
this halo mass range. 

In order to assign a probability to the events shown in Figure 1, we use a Monte-Carlo 
procedure by which we derive R{r) for each of 10^ random halo sub-sample pairs of the 
same size as those of our main analysis (and for each mass range and redshift bin). We 
then ask how many times is R{r) consistently less than 1, for all r's (or for 20 Mpc 
in the case of 2; = 1 and the higher mass bins), to find that in all cases the probability is 
< 0.0009. If we now ask a more restrictive but more accurate question, ie., how many times 
would R{r) be systematically lower at the level observed between the substructured and 
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non-substructured halo sub-samples, then the answer < 10~^. We therefore conclude that 
haloes with substructure are significantly more clustered than haloes without substructure, 
locally and in high redshifts, in agreement with: 

(a) the observational results based on clusters of galaxies (Plionis & Basilakos 2002), 

(b) the DM halo occupation number-clustering correlation (Wechsler et al. 2006; Gao & 
White 2007), and 

(c) the weak signal found by Wetzel et al. (2006) that haloes which have undergone a 
recent major merger (ie., have significant substructure according to our nomenclature) show 
a slightly enhanced clustering. 

To test whether the dependence of the clustering length, tq, on halo substructure is 
due only to those haloes which experience "major" mergers (as suggested by Wetzel et 
al. 2006) we have divided our haloes, of each mass range, in 4 equal-number sub-samples 
based on their sorted A-deviation values (corresponding to different substructure amplitudes) 
and computed their individual two-point correlation function. In Figure 2 we present the 
corresponding correlation lengths as a function of A for all halo sub-samples (the indicated 
A-deviation value corresponds to its median value in each of the 4 ranges) . It is evident that 
the correlation length is growing smoothly and monotonically with increasing A-deviation 
value, a fact which argues against the tq — A relation being only due to "major" mergers. 
A similar result has been found also for real clusters of galaxies (Plionis & Basilakos 2002 - 
their figure 3). 

Now we ask what could the reason be for such an effect? Could it be that haloes having 
substructure (ie., being dynamically active) are located in high density environments (see 
also Ragone-Figueroa & Plionis 2007) and thus at regions of a relatively early halo formation 
time, while haloes with no substructure are located in low-density regions, ie., at regions of 
a later halo formation time (eg. Sheth & Tormen 2004; Zhu et al. 2006; Marker et al. 2006). 
The higher clustering of substructured haloes indeed points in such a direction, although one 
would have expected that earlier formed haloes would typically be dynamically more relaxed 
than equal-mass later formed ones and would have had less evident substructure features. 
However, the higher merging and interaction rate, expected in high-density regions, could 
disrupt the relatively relaxed halo dynamical state. 

To attempt to answer the question posed, we present in Figure 3 the ratio of the nor- 
mahzed SM/M frequency distributions of the substructured and non-substructured haloes 
(for the two extreme mass ranges, just for economy of space). If both type of haloes traced 
similar overdensities, this ratio should have been statistically equivalent to one. However, it 
is evident that it is significantly higher (or lower) than 1 at larger (or smaller) overdensities. 
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for all halo mass ranges and redsliifts. As we had anticipated, substructured haloes are 
typically found in higher- density regions. 



4. Conclusions 

In this letter we have investigated the relation between DM halo clustering, halo dynam- 
ical state and halo large-scale environment. We identified haloes with a FOF algorithm in a 
dark matter only ACDM simulation and used haloes with M > 10^^ h"^ Mq, identified at 
z = and z = 1. The halo dynamical state was determined by measuring the amount of halo 
substructure using an observationally driven approach. We then calculated the two-point 
correlation function for haloes with high and low levels of substructure, finding that the 
former haloes are significantly more clustered than the latter, while there is also a smooth 
increase of the halo correlation length with increasing amplitude of the halo substructure 
index. 

Finally, we find a highly significant signal by which haloes with high levels of substruc- 
ture are located typically in higher density regions with respect to haloes with low levels of 
substructure and this could be an explanation of our previous results. The higher clustering 
of haloes found in high-density large-scale environments should be expected if haloes collapse 
earlier in such regions and thus have more time to evolve and cluster, while their higher- 
levels of substructure should be probably attributed to the higher rate of halo interactions 
and merging which is present in such high-density regions. 
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by the Consejo de Investigaciones Cientificas y Tecnicas de la Republica Argentina (CON- 
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Fig. 1. — The ratio of the 2-point correlation functions of haloes with and without substruc- 
ture for both redshifts epochs. The continuous lines corresponds to z = 0, while the dashed 
lines to 2; = 1. The errors shown are the propagated quasi- Poissonian ^(r) uncertainties. 
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Fig. 2. — Correlation length, function of halo A-dcviation for the three halo mass 

ranges and for two redshift snapshots {left panel: 2; = and right panel: z — 1). The different 
halo mass ranges are indicated by the different symbols. 
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Fig. 3. — The ratio of the normahzed 5M/M frequency distributions of substructured and 
non-substructured haloes. The 5M/M fluctuations have been evaluated on a grid with a 5 
h-^ Mpc cell -size. Results are shown for two mass ranges and two redshifts. Error-bars are 
propagated individual Poissonian uncertainties. 



